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In nucleus-nucleus collisions at RHIC and LHC, partons produced at high transverse 
momentum can undergo multiple scattering within the collision region prior to fragment- 
ing into hadrons. We have studied the resulting medium-modified fragmentation function 
based on a calculation of the BDMPS-Z medium-induced gluon radiation for a dense, 
expanding medium of small finite extension. Here we explain how the BDMPS transport 
coefficient q which measures the energy density attained in the collision, can be extracted 
from the observed modification of high-p^ hadroproduction. We also comment on the 
significant remaining uncertainties in extracting q from data. 

To leading order in perturbative QCD, high-p^ hadroproduction in proton-proton col- 
lisions is described by the factorization formula 



dpfdy 



Here //^, are the proton parton distribution functions (PDF), dcr'^^~^^/dt is the partonic 
cross section and Dk^h{z, Q^) describes the fragmentation of a parton k into the hadron 
h carrying a fraction z of the momentum. Eq. (|l]) leads to a fair description of the shape 
of high-px hadronic spectra while the normalisation has to be adjusted by an energy- 
dependent i^'-factor (see e.g. [0]). Also to NLO, the disagreement between theory and 
experiment lies essentially in an albeit reduced normalisation factor f^. 

For nucleus-nucleus collisions, additional corrections arise due to parton density effects. 
At sufficiently high pj_ where the factorised form (|l|) is expected to hold, two corrections 
are known: i) Density effects in the initial state can be parametrised by nuclear PDFs 
for which global fits are available [Q]. At RHIC energies, however, this is only a small 
correction |I| [^]. ii) Density effects in the final state affect the fragmentation of the high 
p± parton produced in the hard collision. Several studies indicate that the dominant final 
state correction is due to medium- induced gluon radiation and the ensuing parton energy 
loss. If the high p± parton looses with probability Pe{^) a fraction e of its energy while 
propagating through the medium, then its medium-modified fragmentation function can 
be written as 
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We have calculated the quenching weights Pe{^) for quarks and gluons starting from 
the BDMPS-Z P,p. [T0|JTl[1 gluon radiation spectrum dl/du generalised to the case of an 
expanding medium of small finite size [0]. For the case of a static medium, the resulting 
finite size corrections to the BDMPS-expression of dl/duj are shown in Fig. 1. They 
become important for gluon energies much smaller than the characteristic gluon frequency 
Uc = q L. The latter depends on the in-medium path-length L and the BDMPS transport 
coefficient q which measures the amount of transverse momentum squared transferred 
from the medium to the hard parton per unit path-length. For a static medium, the 
radiation spectrum depends only on Uc and the dimensionless parameter combination 
R = \qL^. As seen in Fig. 1, the BDMPS case is the limit L ^ oo of our calculation 
keeping ojc finite. 
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Figure 1. (Up) Gluon radiation spectrum, 
Lodl/du, for several values of R and compar- 
ison with the BDMPS result. 

Figure 2. (Right) Continuous (upper fig- 
ure) and discrete (lower) contributions to the 
quenching weight Pe{^) of a fast quark for 
different values of R. 




Quenching weights are evaluated in the independent gluon emission approximation |]T2 
and show both a discrete and a continuous contribution (see Fig. 2), 
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Here, -P(e) is a generalised probability of norm 1. For a medium of finite length L, 
there is a finite probability po that the fragmentation is not disturbed by the medium. 
This discrete contribution decreases with increasing in medium path-length or increasing 
density of the medium. 

The transport coefficient is proportional to the density of scattering centres, g(r) oc 
n(r). The dynamical expansion of the collision region can be parametrised by a power 
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law decrease in proper time, ?7,(r) = no (^^) , where a = for static, a = 1 for Bjorken 
expansion. Fig. 3 shows the gluon radiation spectrum for different values of a with 
parameters Uc and R written in terms of the linearly weighted transport coefficient 
2 r-ro+L 
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dr (r - To) g(r) . 
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After rescaling according to (^, the radiation spectra dl /duo calculated for different values 
a coincide (see Fig. 3). This establishes a dynamical scaling law [0] which simplifies cal- 
culations considerably since results for a dynamically expanding medium can be obtained 
by scaling the result obtained for a static medium. It extends earlier findings of dynamical 
scaling of the average energy loss ||I0|J6|JT^ to the tu-differential radiation spectrum 0. 
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Figure 3. The gluon radiation spectrum, 
udl/duj, plotted function of lUc = qL 

and R = \(lL^- Different curves correspond 



to different expansion parameters a 
1.0 and 1.5. 
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Figure 4. The LO KKP [14] fragmenta- 
tion function u ^ it for no medium and 
the medium-modified fragmentation func- 
tions for different gluon rapidity densities 
in a medium of length L = 7 fm. 



Based on the quenching weights (H) for static and expanding media, we have calculated 
medium- modified fragmentation functions according to (Q), using the KKP parametrisa- 
tion [|I4| for parton fragmentation in the vacuum (see Fig. The virtuality Q of the 



parent parton is identified with its initial transverse energy. For a medium showing one- 
dimensional Bjorken expansion (a = 1), the transport coefficient q is related to the gluon 
rapidity density as 

1- , L2 dN9 

The calculation presented in Fig. ^ allows to determine medium modifications to the 
transverse momentum spectrum (|lD at sufficiently high p±. While this last step is not 
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yet performed, one can estimate the effect based on the observation |jT| that the hard 
partonic cross section in (|1]) essentially weighs the fragmentation function with the power 

. Inspecting the corresponding unintegrated moment z^D^jJ^z^ ~ p^) in Fig. ^, one 
estimates that a factor 2 suppression of the hadronic pion spectrum for ~ 6-7 GeV 
corresponds to an initial gluon density of 500-1000. For an in-medium path-length of 
L = 7 fm, this amounts to a line-averag ed transport coefficient q = (300 MeV)^ - (600 
MeV)^/fm which corresponds to (pi)/i^med of several GeV^/fm within the first two fm/c. 

For a reliable determination of q from data, the evaluation of (|l]) is not sufficient. One 
also has to discriminate the medium modification discussed here from other corrections 
to (|l]), such as "soft" hadronic contributions at low p± and higher order perturbative con- 
tributions at high p±. Much of this can be done by studying the dependence of spectra 
on the in-medium path-length e.g. via the centrality dependence of high-p^ hadroproduc- 
tion and monojet/dijet production. Another discriminatory tool results from the fact that 
medium-induced gluon radiation is a medium-enhanced power correction which shows a 
power law decrease with increasing p±. This is in contrast e.g. to perturbative NLO 
corrections which mainly affect the normalisation but not the shape of spectra. This will 
also provide an important additional cross check at LHC: according to (^), if the rapidity 
density increases by a factor 5 from RHIC to LHC, and if the suppression factor 2 for 
RHIC 7r° spectra at p± ~ 6-7 GeV can be attributed solely to medium-induced gluon 
radiation, then a factor 2 suppression of the corresponding spectrum at LHC is expected 
in the range p± ~ 40-50 GeV. At even higher p± still accessible at LHC, this suppression 
would decrease according to a power law, thus leading to a modified p^-shape of hadronic 
spectra. 
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